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Abstract 

We consider the problem of inserting a stiff chain into a colloidal sus- 
pension of particles that interact with it through excluded volume forces. 
The free energy of insertion is associated with the work of creating a cav- 
ity devoid of colloid and sufficiently large to accommodate the chain. The 
corresponding work per unit length is the force that resists the entry of 
the chain into the colloidal suspension. In the case of a hard sphere fluid, 
this work can be calculated straightforwardly within the scaled particle 
theory; for solutions of flexible polymers, on the other hand, we employ 
simple scaling arguments. The forces computed in these ways are shown, 
for nanometer chain and colloid diameters, to be of the order of tens of 
pN for solution volume fractions of a few tenths. These magnitudes are 
argued to be important for biophysical processes such as the ejection of 
DNA from viral capsids into the cell cytoplasm. 



1 Introduction 

DNA is known to be packaged at very high concentrations in most viruses pQ. 
In the particular case of bacterial viruses, or bacteriophages ("phage"), the 
high free energy density of double-stranded DNA bent and crowded inside the 
viral capsid is thought to be responsible for the initial injection step of the 
viral genome into its host cell. In these instances the virus binds to its specific 
receptor on the outer cell membrane, its capsid is opened, and the stored energy 
of confinement drives the ejection process. As the ejection proceeds the driving 
force decreases because of the progressive relief of stress inside the capsid. Later 
steps in the translocation of the full viral genome from the capsid into the cell 
vary from one type of phage to another. For example, it has been shown in the 
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case of T7 bacteriophage that the translocation is coupled to the transcription 
of viral DNA by the host cell machinery, which exerts a force pulling the DNA 
inside the cell In other cases, such as T5, the translocation can occur without 
the intervention of an active process, e.g., via interaction with the DNA-binding 
proteins present at high concentrations in the cell cytoplasm [3]. 

Over the last couple of decades there have been several attempts to estimate 
theoretically the energy cost of confining DNA inside a viral capsid ^|SJ|2]. 
Recent experiments [7| and simulations and phenomenological theory [51I9HTU] 
on the packaging of viral DNA indicate that the stress associated with fully 
loaded capsids corresponds to ejection forces on the order of tens of pN and to 
internal pressures of several tens of atm. Osmotic pressures of this magnitude 
are commonly realized in polymer solutions at moderately high concentration, 
thereby raising the following simple question [IT: is it possible to stop the 
DNA ejection from a viral capsid by "opposing" the internal pressure with 
that exerted by a colloidal suspension? This question can be formulated more 
precisely in terms of the associated forces: can a colloidal suspension (hard 
spheres or flexible polymers) produce a force opposing that driving the DNA 
out of the capsid? This problem is also relevant for in vivo phenomena, since the 
cell cytoplasm, containing high concentrations of various components (mainly 
proteins), has been modeled as a colloidal suspension to evaluate the effect of 
the "macromolecular crowding" on various biochemical reactions |12) . 

Our aim in this work is to formulate a theory for the force resisting chain 
insertion into a colloidal suspension. This force is evaluated by calculating 
the free energy of insertion of the chain in the suspension, more explicitly the 
work required to "grow" the particle to its final size and shape in the solu- 
tion PS]- This concept has long played a fundamental role in the evolution 
of modern theories of simple liquids [II]. For the case of hard sphere fluids, 
one of the most physically appealing and natural tools for calculating this work 
of insertion is the Scaled Particle Theory (SPT) |15|. The force resisting the 
injection is then given by the derivative of this work with respect to inserted 
chain length. If the colloidal suspension is a solution of flexible polymers, the 
work of insertion is evaluated naturally through scaling laws. As will be shown 
in this paper, the force resisting chain injection in both of these instances can 
indeed be comparable to the force driving DNA out of viral capsids |51l^lll0|. 
This result accounts for recent experiments in which DNA ejection from bac- 
teriophage (A) is shown to be suppressed by a sufficiently high concentration 
of flexible polymers (polyethyleneglycol - PEG) in solution [lfij . In these mea- 
surements an enzyme is also present which digests the ejected DNA into free 
nucleotides, so that one does not have to deal with any of the complications of 
PEG-induced DNA condensation ^7|, since only short (sub-persistence length), 
isolated, pieces of DNA are involved. This also validates our modeling of the 
ejected DNA as a stiff chain, in calculating the force resisting its entry into 
the colloidal suspension. More systematic analyses of these experiments will be 
provided elsewhere. Here we focus on the relationship between resisting force 
and osmotic pressure, comparing and contrasting the situations of hard sphere 
versus flexible polymer solutions. 
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Stiff chains - rods - in both hard sphere suspensions and flexible polymer 
solutions have been treated previously, using SPT and scaling theory, respec- 
tively. But the aims have invariably been the calculation of phase behavior - 
say, demixing via depletion interactions in rod/sphere mixtures or chain 
condensation via effective attractions in solutions of flexible polymers • Our 
aim, on the other hand, is to describe directly the forces that arise in these 
colloidal systems, acting directly to resist the entry of chains. The paper is 
organized as follows. In the next section the general principle of the SPT is 
outlined briefly and results from earlier work are quoted for the calculation of 
the work of insertion of a spherical cavity in a hard sphere fluid. In section|21wc 
generalize this approach to the case of a cavity whose shape is spherocylindrical, 
calculating the force resisting entry of the chain from the insertion energy per 
unit length. The work of inserting a cylinder into a flexible polymer solution 
is then estimated from simple scaling arguments in section 0] again generaliz- 
ing from the spherical solute case. The forces derived in both instances - hard 
sphere and flexible polymer solutions - are discussed in the final section, as a 
function of chain diameter, colloid size, and concentration. Numerical estimates 
are made for the case of DNA and PEG solutions, before we conclude with 
some remarks about connections to other biophysical phenomena involving the 
translocation of stiff chains into concentrated colloidal suspensions. 

2 Work of insertion of a spherical cavity in a 
hard sphere fluid 

We present in this section the SPT calculation of the work of insertion of a 
spherical cavity of arbitrary radius r in a hard sphere fluid of concentration cp, 
the diameter of the spheres being a. For the sake of simplicity, we choose kT 
as the energy unit throughout this paper. The cavity is defined by a spherical 
volume free of any centers of hard spheres. For a small enough cavity, 
r < a/2, the work of insertion is related exactly to the free volume fraction by 



since there is only room for zero or one sphere. For bigger cavities, more spheres 
can be accomodated, so that exact determination of W(r) would require the 
knowledge of n-body correlation functions 15 . 

One of the key achievements of the SPT is the providing of a physically mo- 
tivated approximation for the work of insertion without determining the whole 
structure of the hard sphere fluid. Within the SPT it can be shown that the form 
of the work of insertion of a spherical cavity is restricted by a certain number of 
exact conditions. By choosing a simple form for the work of insertion satisfying 
those exact conditions, one expects to have an accurate approximation for this 
work. Since we want to calculate the work of insertion for an arbitrary value of 
cavity radius, it is judicious to choose the asymptotic value for r — > oo given by 




(1) 
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the sum of a pressure-volume work and a surface tension contribution: 



W 00 (r)=pV + a[l-—)S (2) 



r 



The quantities p, a, S are respectively: the pressure of the hard sphere fluid; 
its surface tension in the limit of flat surface; and the Tolman length associated 
with the curvature of the surface [23j- V — 4irr 3 /3 and S = 4irr 2 are the volume 
and surface of the cavity. As shown below, by using three conditions derived 
exactly through the SPT, one can determine the hard sphere fluid properties. 

The first two conditions are the continuity of the work's first two derivatives 
at r = a/2. Therefore, by matching the macroscopic values of these quantities - 
see Eq(3 with their exact microscopic values for r < a/2 - see EqQ] we obtain 

V(l + y) 



2-Ka 2 {l-yY 
and 

y = 7ra 3 cp/6 being the volume fraction of spheres. Note that the surface tension 
of a hard sphere fluid is negative, reflecting the absence of attractive interactions 
holding the particles together as in a "real", cohesive, fluid. The pressure is 
found by another exact relation, 

^ = l + 4yg(a), (5) 
cp 

where g(r) is the pair correlation function, proportional to the probability of 
finding the center of a hard sphere at a distance r from the origin, provided 
that there is already another particle at the origin |14j . 17(a) can in turn be 
expressed, without approximation, in terms of the first derivative of the work 
of insertion: 

In this way one obtains the SPT equation of state for the hard sphere fluid 

= _6_ y + y 2 + y 3 , , 

P tto 3 (1-y) 3 {) 

The second and third virial coefficients derived from the expansion of Eq0 for 
low y are known to be exact, while up to the fifth order the discrepancy is less 
than 5 %. Therefore as long as the volume fraction of particles is not too high, 
the SPT equation of state provides an accurate description of the hard sphere 
fluid. Eq. is also known to be identical to the pressure equation of state 
obtained via the Percus-Yevick closure 21 , 2211231. 
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3 Work of insertion of a spherocylinder in a 
hard sphere fluid 

We generalize in this section the SPT formalism to calculate the work of in- 
sertion of a spherocylindrical cavity in a hard sphere fluid. The spherocylinder 
considered here has a total length La + (1 + D)a and a diameter (1 + D)a. 
When L, D = 0, the cavity is just a sphere of diameter a, so that the exact 
form of the work of insertion is still given by Eq^ For bigger cavities, we can 
compute it approximately by using, as in the original SPT approach, exact con- 
ditions restricting its form. A quite natural choice is the work of insertion of a 
macroscopic spherocylinder in a hard sphere fluid 

W(D,L)=pV + a 1 (l-^St + ai (l - |Q S 2 (8) 

We explicitly choose different coefficients associated with the spherical (indices 
1) and cylindrical (indices 2) part of the cavity, in order to start with the most 
general form of the macroscopic work. As shown below, these quantities turn 
out to be equal in the present calculation. i?i and i?2 are respectively the 
curvature radii of the spherical and cylindrical parts of the cavity. This form of 
the work of insertion again corresponds to a pressure- volume work and a surface 
contribution. Note also that this form involve an expansion of the work in powers 
of D and L, since V = to 3 ((1 + Df + + D) 2 ) /6, S x = 7ra 2 (l + D) 2 and 
S*2 = ita 2 L{\ + D) [211118) . The coefficients are found by matching the first two 
derivatives of the work at D,L = 0. The results read 

3y(2 + j/) pa 

ai=(T2= 2na 2 (l-y) 2 ~Y (9) 

and 

Sl = 52 = i ^Tlfeir ( 10 ) 

p7ra (i-y) 2 

Finally, the pressure involved in the preceding equations is found by using Eqs|5] 
and El leading again to expressions for p, a and 5 given by Eqs0|Hland2] The 
fact that the coefficients are simply the pressure, surface tension and Tolman 
length of a pure hard sphere fluid indicates that, for such a simple geometry, 
the shape of the cavity does not really matter for the work of insertion, and the 
macroscopic form Eq. [5]obtains as before. The work of insertion Eq|S| together 
with EqEElandU is the same as the one derived previously for phase behavior 
studies of rod/sphere mixtures 18 ,23). The work of insertion in the present 
study will be used in section to evaluate the force resisting the insertion of a 
spherocylinder into a hard sphere fluid of given volume fraction. 
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4 Work of insertion of a cylinder in a flexible 
polymer solution 

Unlike in the hard sphere fluid case, no exact conditions restricting the form of 
the actual work of insertion of a cavity of arbitrary shape or size in a polymer 
solution has been derived yet. Moreover it is not straightforward to try to map 
the SPT principle to the polymer case because the cavity created around one 
chain excluding other chains is itself a statistical object, i.e. depends on the 
conformation of the polymer. Nevertheless, the work of insertion of a cavity 
of spherical or cylindrical shape can be calculated using either pure scaling 
arguments |25H26llT§] or renormalization group arguments |27II28| . i.e. specific 
polymer physics tools. We briefly present here scaling arguments that allow one 
to approximate the work of insertion of a cylindrical cavity in a polymer solution 
by superimposing asymptotic behaviors of the work. This means implicitly that 
all numerical prefactors will be set equal to unity. 

We consider first the interaction of a cylinder of length L and diameter D 
with a flexible chain of radius of gyration R g in the limit of thin cylinder R g » 
D. The monomer size and the degree of polymerization are, respectively, b and 
M. The interaction energy is proportional to the number of contacts between 
the cylinder and polymer. If the two objects occupy the same average volume 
i? 3 , then their interaction energy within a simple mean-field approximation is 
given by 

~ R? <\>cyl<\>coil (11) 

where the quantities <p cy i ~ ^p- and 4> co a ~ ^pj- are respectively the volume 
fraction of the cylinder and the chain. Therefore we expect this interaction 
energy to be linear in both the length of the cylinder and the degree of poly- 
merization of the coil. The mean-field argument used here is justified since 
the interaction energy scales with the common size of the interacting objects 

7-)2. 4/3 

R as F int ~ kT R ° /3 , and therefore decreases as the sizes of the objects are 
increased j^Sl- Similar arguments can be used to show for example that the 
contribution of binary contacts in a single chain is irrelevant if the dimension of 
the embedding space d exceeds 4, while in the opposite case simple mean field 
approximations fail 29 . 

It is straigthforward to show that regardless of whether the cylinder is longer 
or smaller than the radius of gyration of the chain, the interaction energy must 
still vary linearly both with the length of the cylinder and the degree of polymer- 
ization of the chain. With these requirements, a cross virial coefficient between 
cylinder and polymer of general scaling form v cy i-ch ~ L a D°R g ' is uniquely 
determined and becomes 

Vcyl-ch ~ LD^R 5 / 3 (12) 

In a dilute solution of polymers of monomer concentration c, the work of inser- 
tion is then given by 

W^LD^R 5 / 3 ^ (13) 
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The previous arguments can easily be generalized to the case of semi-dilute 
solution of polymers where chains overlap. The overlap concentration is given 
by c*b 3 — y* ~ Af~ 4//5 [2H]. In this limit, the polymer solution is characterized 
by the average mesh size of the temporary network formed by the chains, or 
equivalently the correlation length, that scales like £ ~ 6(c6 3 )~ 3 / 4 A chain 
subunit of size £ is called a semi-dilute blob, so that the semi-dilute solution 
can be envisioned as a compact packing of such blobs, with the statistics of 
monomers within a blob given by those of an isolated chain: if g is the number 
of monomers per blob, then £ ~ g 3 / b b. The cross virial coefficient between the 
cylinder and the blobs is given by v cy i-u ~ LD 1 ^ 3 ^ 5 ^ 3 , so that the work of 
insertion in a semi-dilute solution reads 

W~LDV 3 e /3 ~~^ (14) 

As noted by de Vries, the result for the work of insertion in the dilute and 
in the semi-dilute case are the same: 

LD 1 ' 3 , . 

w ~-^ v ( 15 ) 

with y = cb 3 |19|. This work is independent of the chain degree of polymeriza- 
tion, and is linear in the volume fraction of monomers. This follows from the 
fact that the insertion of a cylinder is dominated by the depletion layer around 
it of order D and is a local process, as long as its diameter D is smaller than 
the correlation length of the polymer solution (R g in dilute solution and £ in 
semi-dilute solution). 

When R g , £ << D on the other hand, the depletion layer around the cylinder 
is of order the correlation length of the polymer solution, so that the work of 
insertion is dominated by the pressure- volume work and surface tension contri- 
butions. As in the case of a large cavity in the hard sphere fluid, this corresponds 
to the macroscopic value of the work. Therefore this work is written 

f ^jLD 2 + a d LD if c « c* 
W~{ ^ LD* + a sd LD if c»c* ^ 

The surface tensions scale like a d ~ jjR g and u sd ~ 28 . Therefore the 
different asymptotic behaviors of the work of insertion of a cylinder in a polymer 
solution can be summarized as follows. 



if R g < D, 



if R g > D, 




„- w— ■ %&DLfoTy<y* 

W rod -\ _ a/ALD , , ym nL fQry>y * (17) 



W *= \ y "Hf) 1/3 fOT 2/ < 2/i 

rod { y 9 ^ + y 6 ^ior y>yi 
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where the volume fraction yi ~ (j)) 4 ^ is defined by the condition £ = D. One 
can check that the condition y* < yi is equivalent to D < 

Since the rod (DNA) and flexible polymer (PEG) used in the experiments 
of reference ^HI satisfy D ~ corresponding to the crossover of the previous 
expressions of the work of insertion, we can try to approximate the actual work 
of insertion by simply superimposing the asymptotic behaviors, following de 
Vries . Therefore the general form of the work of insertion of a cylinder in a 
polymer solution is given by 

W = viy + v 2 y m + v 3 y 9/4 (19) 

with the coefficients 




•* ~ -yf (21) 

, 3 ~ ±f (22) 

These results are identical to those of de Vries, except for our having taken into 
account explicitly the surface tension of the polymer solution in order to be 
consistent with the SPT approach. 



5 Numerical estimates 

We calculated in the last two sections the work of insertion of a cylinder (or 
a spherocylinder) in a colloidal suspension in both the cases of hard spheres 
and polymers. The classical interpretation of this quantity is the energetic cost 
associated with the growth of the particle to its final size and shape in the bulk 
suspension. This work is therefore associated with a characteristic force that 
resists the growth of the particle. In the case of the ejection of DNA from a viral 
capsid into a colloidal suspension, the DNA is "inserted" into the suspension 
along its length. Therefore the force associated with the work of insertion is 
simply given by the derivative of the work with respect to the length of the chain. 
In the experiments done by Evilevitch et ai, the DNA is digested by the enzyme 
DNase I as it is inserted from the capsid into the colloidal suspension, suggesting 
that we approximate the ejected DNA as a straight cylinder of diameter (Idna- 
This allows us to neglect the effect of DNA condensation and aggregation. In the 
case of a hard sphere fluid, the force is calculated by using the SPT expression 



for the work of insertion of a cylinder: 



dW(D, L) 
d{La) 



The prefactor 4 gives the force in pN, provided that both the diameters of the 
cylinder d^NA and of the spheres a is expressed in nm. Note that a sphero- 
cylindcr of diameter doNA in & hard sphere fluid creates a cavity of diameter 
dDNA + a excluding the center of any hard spheres. The first term in Ea l23l 
is the pressure contribution while the second is the surface contribution. The 
magnitudes of those two contributions are shown in figure ^ for a particular set 
of diameters. It is clear from this figure that a simple argument where the work 
of insertion is approximated by a pressure- volume work term overestimates the 
force, but this effect is not dramatic. The influence of the sphere size is shown 
in figure |3 The force gets smaller as the sphere size is increased at fixed volume 
fraction - fewer particles contribute to the osmotic pressure. If we think of 
various proteins found in vivo behaving purely as hard spheres |12| . the range 
of their sizes from 1 to lOnm shows that for moderate colloid volume fraction 
(typically 20-30% ^2) the magnitude of the force resisting insertion is of the 
same order of magnitude as the force driving DNA out of a viral capsid, mea- 
sured and estimated [HHH1E31 recently to be tens of pN in the first steps of 
DNA ejection. 

Similarly in the case of polymer solution, the force resisting ejection is given 

by 

/ = ^y + 4y 6/4 + ^y 9/4 (24) 

with v[ — Vi/L, the i/j's being given by Eas l2QI22l Forces / vs volume fraction y 
is shown in figure|31for different choices of monomer sizes. The value b = QAnm 
approximates the statistical size of ethylene glycol monomer 30 . As mentioned 
earlier, the work of insertion of DNA in a solution of PEG has already been 
calculated by de Vries in his effort to address the polymer-salt-induced DNA 
condensation problem; his fit of the numerical prefactors to experimental data 
gives an f(y) close to ours for b = OAnm. Once again, as in the hard sphere 
case, the force is of order tens of pN, comparable to the force ejecting DNA 
from phage capsids (71819X101 . Preliminary experimental results show that DNA 
ejection from A-phage can indeed be inhibited by a PEG solution at moderate 
concentrations, the force pushing out the DNA being balanced by the force 
resisting ejection JSj. At high enough polymer concentration, no ejected DNA 
is detected. The smallest concentration required for the complete suppression 
of DNA ejection, up to experimental limits of detection, implies a resisting force 
between 20 and 30 pN. This corresponds to the right order of magnitude of the 
force driving the DNA out of the capsids when it is fully loaded, as estimated 



3 1 



dp? 



(y + . 

2a(l-y) 3 



V) 



-4^(1+^(1+3/)) P N 
2a{l-yf 



(23) 
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by a careful balance of bending energy and intermolecular repulsions inside the 

capsid miaun]. 

6 Concluding remarks 

We calculated in this study the work of insertion of a cylinder in a hard sphere 
fluid and in a polymer solution of arbitrary concentration. In the context of 
DNA ejection from a viral capsid, this work can be interpreted as the work 
done against the constant force resisting insertion of the DNA in the suspension. 
This allows us to evaluate this force by differentiating the work of insertion with 
respect to chain length. 

In the case of hard spheres, the work of insertion is computed using the 
SPT, which is known to provide accurate thermodynamical properties of the 
hard sphere fluid provided that the volume fraction is not close to the close- 
packing limit. It is shown that the order of magnitude of the work of insertion 
can be approximated simply through a pressure-volume work term. The force 
computed in this way is of the order of several pN for sizes of hard spheres 
ranging from lnm to lOnm. This range of forces compares favorably to the 
order of magnitude estimated for the ejection force of DNA from bacteriophage. 

In the case of flexible polymer solution, the work of insertion is calculated 
using scaling arguments; the force estimated in this way depends on omitted 
numerical prefactors which we drop, an omission inherent in the scaling meth- 
ods. The order of magnitude of the work of insertion can still be evaluated by 
setting those prefactors equal to unity. We estimate the resisting force exerted 
by a PEG solution to be of order tens of pN for moderate PEG concentrations. 
This force is again comparable to the force driving DNA out of capsids |51l51ll(J| . 
a result that has been qualitatively confirmed by recent experiments in which it 
is shown that the amount of ejected DNA decreases with PEG concentration. A 
more thorough analysis of those experiments will be presented elsewhere using 
the present theory . 

The spirit of the SPT calculation for the work of insertion in the case of 
hard spheres and that of the scaling calculation in the case of polymers are very 
similar. In both cases, the functional form of the work of insertion and its prac- 
tical evaluation are guided by mathematical and physical requirements. The 
mathematical requirement in the case of the SPT is the choice of an analytical 
function to approximate the work, motivated by the appearance of disconti- 
nuities only in high order derivatives. In the case of polymers, scaling laws 
are required because of the critical properties of polymer solutions in the in- 
finite chain length limit. In both cases, the work of insertion matches certain 
asymptotic and well-understood behaviors that are physically motivated. 

The present calculation provides also a concrete realization of the concept 
of work of insertion, and of the force associated with it; this concept underlies 
many modern theories of fluids, but now finds a direct manifestation in the 
experiments where DNA is "inserted" into a colloidal suspension from a viral 
capsid. 
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Figure Captions 



• Figure 1: Different contributions to the force resisting ejection in pN, as 
a function of volume fraction. The full line is the total force, the long 
dashed line is the pressure contribution and the short dashed line is the 
negative surface contribution. cLdna = 2nm , a = lnm. 

• Figure 2: Force resisting ejection in pN, as a function of volume frac- 
tion, for different colloid diameters. From top to bottom a = Q.5nm, a = 
lnm, a = 1.5nm, a — bum. djjNA = 2nm. 

• Figure 3: Force resisting ejection in pN in a PEG solution, as a function 
of PEG volume fraction for different monomer diameters. From top to 
bottom b = 0.2nm, b = OAnm, b — 0.6nm. doNA — 2nm. 
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Figure 1: Different contributions to the force resisting ejection in pN, as a 
function of volume fraction. The full line is the total force, the long dashed line 
is the pressure contribution and the short dashed line is the negative surface 
contribution. cLdna = 2nm , a — lnm 
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Figure 2: Force resisting ejection in pN, as a function of volume fraction, for 
different colloid diameters. From top to bottom a = 0.5nm,a — lnm,a = 
1.5nm,a = 5nm. djjNA = 2nm 
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Figure 3: Force resisting ejection in pN in a PEG solution, as a function of 
PEG volume fraction for different monomer diameters. From top to bottom 
b = 0.2nm, b = OAnm, b — 0.6nm. cLdna = ^nra. 
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